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Half-quantum vortices (HQVs) can exist in a superconductor or superfluid with an exact or ap- 
proximate U(l)xU(l) symmetry, for instance in spinor condensates, 3 He-A, Sr2RuC>4, and possibly 
cuprate superconductors with stripe order. In this paper, we show that a lattice of HQVs can be 
stabilized at finite temperature even when it does not have lower energy than the lattice of full 
vortices at T = since there is a gain in configurational entropy when a full vortex fractionalizes 
into a pair of HQVs. Specifically, the lattice of HQVs has an optical branch of phonon modes absent 
in the lattice of full vortices. Moreover, the HQV lattice at T > can have a different structure 
than the HQV lattice at T = 0. 



I. INTRODUCTION 



Superconductors and superfluids which require two 
phase variables to describe their condensates can be said 
to possess a U(l)x U(l) symmetry. One example, which 
has been known for a long time, is a thin film of the A- 
phase of 3 He superfluid 1 2 in a perpendicular magnetic 
field, where the spin-triplet Cooper pairing allows for 
the possibility of the Cooper pair spin rotating in the 
plane perpendicular to the field. More recently, Sr2Ru04 
has been identified as an electronic analogue of 3 He-.ApEI 
and as such can allow for a U(l)xU(l) symmetrjEK^. 
A U(l)xU(l) symmetry can also arise in Bose conden- 
sates of atoms with nonzero integer spins, commonly 
known as 'spinor condensates EHIS] Lastly, "pair-density- 
wave" (PDW) order - long known as th e Fulde-Ferrell- 
Larkin-Ovchinnikov (FFLO) phase^EI] - can also al- 
low for a U(l)xU(l) symmetry as the superconducting 
and CDW phases are intertwined there. Recent, it was 
proposecpSHUl t na t experimental studies of the effects of 
the unidirectional spin and charge ordering ('stripe or- 
der') in the cuprate su perco nductors La2-iBa K Cu04 and 
Lai.e-zNdo^Sr^CuOi 121221 shows evidence for the PDW 
order. Another candidate system for a PDW order is 
CeCoIns under strong magnetic fiel d 21 * 22 *. Possible phase 
transit ions i n various systems with PDW order has been 
studiecP^HUl. 

A condensate with a U(l)xU(l) symmetry allows for 
the existence of a half-quantum vortex (HQV) with only 
7T, rather than 2ir, winding in the overall phase of the 
superfluid, 6. For a charged superfluid, i.e. a supercon- 
ductor, a HQV binds flux of h/Ae, which is half the flux 
of the conventional Abrisokov vortex, $o = h/2e. This is 
a consequence of the existence of an extra phase variable 
a, where at long distances, the condensate consists of 
two components with phases (9 ± a)/2, regardless of the 
physics that gives rise to a. The single- valuedness of the 
condensate wave function around a vortex is thus main- 
tained even with A8 — ±ir as long as it is also accompa- 
nied by a ±7r winding in a, thus making a HQV topolog- 
ically stable. Conversely, observation of HQVs indicates 
the existence of the additional phase degree of freedom 
in the condensate, and thus serves as an indicator of the 



structure of the order parameter. The HQVs in chiral 
triplet superconductors have attracted particular inter- 
est in recent years stemming from the possibility that 
they will exhibit the simplest non-Abelian statistics^HUl 
More complex non-Abelian statistics of HQVs of spin-3/2 
spinor condensate has also been studied^. 
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FIG. 1. Two possible vortex lattice phase diagrams are possi- 
ble, depending on the ratio of core energies; T is temperature 
and k and p are the phase stiffnesses for a and 9, respectively. 
Note that without any change in the phase stiffness ratio n/p, 
the HQV lattice can be stable at non-zero temperature even 
when it is not stable at T = 0. (n/ p) c is increased if the core 
interaction between HQVs is repulsive, so the upper phase 
diagram is for the case that the vortex core interaction favors 
HQVs over full quantum vortices. 

Although they are topologically stable, this does not 
guarantee that HQVs will ever occur in equilibrium sys- 
tems. Indeed, in superconductors, HQVs have never been 



detected in bulk, due to the energetics issues pointed out 
in Ref[7l To the best of our knowledge, the recent ap- 
parent experimental detection of half-quantum fluxoids 
in mesoscopic Sr2Ru04 sampled^ is, to date, the only 
successful experimental observation of a HQV in a single- 
crystal superconductor. However, in a lattice of HQVs, 
the lattice constant effectively plays the role of the system 
size, so a system close to H C 2 effectively becomes meso- 
scopic, as was shown in Ref. EH The T — structure of 
the HQV lattice obtained for this case is essentially the 
same as that obtained for spinor condensates in Refs. 
ITT1 and H31 It is also analogous to the vortex-antivortex 
lattice configurations proposed in Refs. [53] and EH for a 
two-dimensional (2D) superfluid. 

The main result of our paper is the vortex lattice phase 
diagram shown in Fig. ([T]) that shows the HQV lattice 
favored by non-zero temperature. Studies of the stabil- 
ity of the HQV lattice previously have focused on the 
T = case, i.e. on energy rather than free energy. There 
has not been a thorough investigation of the possibility 
that the HQV lattice may be stabilized by its larger en- 
tropy. While Barnett, Mukerjee and Moore considered 
the possibility of a finite temperature transition from a 
full vortex lattice to a fractionalized vortex lattice^, they 
only treated the case in which this is a result of the tem- 
perature dependence of the ratio of the stiffnesses of the 
two phase variables. Although this is ultimately also an 
entropic effect, it is special to the case in which a is as- 
sociated with spin degrees of freedom with only weakly 
broken SU(2) invariance. In the present work, we will 
show that regardless of any physical origin of a, entropy 
favors the HQV lattice, giving rise to a finite tempera- 
ture phase transition from the full vortex lattice to the 
HQV lattice. In addition, we find the possibility that 
the entropy stabilized HQV lattice can have a different 
structure than the HQV lattice at T = 0. 

The paper is organized as follows: In Section II, af- 
ter comparing the free energy of a full vortex with a 
pair of HQVs, we show that a HQV lattice has higher 
entropy than a full vortex lattice by analyzing lattice 
phonon modes. In Section III, we show the energy com- 
petition between the triangular full vortex lattice and 
HQV lattices with different structures. In Section IV, we 
demonstrate that the HQV lattice at finite temperature 
is always dynamically stable in the thermodynamic limit. 
We conclude with a discussion in Section V. 
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FIG. 2. Three possible vortex lattice structures. The trian- 
gular lattice of full vortices is shown on the top. A full vor- 
tex can fractionalize into a pair of HQVs of opposite types, 
represented by a red dot inside a blue circle. Each type of 
HQV forms a triangular lattice, but depending on the rel- 
ative equilibrium displacements between the two types, the 
two different HQV lattice structures shown in the middle and 
the bottom are found. 



where M* is the Cooper pair mass, A = is the 

London penetration depth, and p and n are the stiffness 
for 9 and a, respectively. In this minimal model, we 
ignore any consideration of crystalline symmetry of the 
superconductor. We can treat the case of a neutral su- 
perfluid by setting e — > (and consequently A — > oo) in 
all results. 



II. HQV ENTROPY 

Results of this paper are based on the minimal model of 
a superconductor in the London limit with a U(l)xU(l) 
symmetry: 



H[6,a] = 
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(1) 



A. Entropy argument 

We first show that a full quantum vortex with 
(A8,Aa) = (2ir, 0), that is the 2-7T phase winding only 
in the overall phase, can always lower its free energy by 
fractionalizing into a pair of HQVs, one with (AO, Aa) = 
(ir,n), the other with (AO, Aa) = (n,—n), because 
this fractionalization will increase the configurational en- 
tropy. We emphasize here that this holds even in the case 
where the pair of HQVs have higher energy at T — 0. As 
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far as T = energy is concerned, in the limit A —> oo 
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where i? is the system size, r is the separation between 
two HQVs, and £ is the core radius, and we assume R ^> 
r 3> £ (in which limit the core energies are negligible). 
The pair of HQVs have higher energy if k > p. However, 
dissociation of a full quantum vortex into a pair of HQVs 
results in an entropy increase 



As(r) = In 
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from the number of the possible HQV pair configurations 
with mean separation r. This implies that even when 
k > p, the free energy / = e — Ts of the pair is lower at 
temperatures Tkt ^ T > T c [ is 
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where Tkt is the Berezinskii-Kosterlitz-Thouless (BKT) 
temperature. 

This observation motivates us to compare the free en- 
ergies of the full quantum vortex lattice and the HQV lat- 
tice in the case in which the latter has higher ground state 
energy. The regime we will examine is H c \ <C H <C H C 2, 
where screening of the charge current (i.e. the finite size 
of A) is not a dominant effect but all vortex cores are still 
well separated. 

We begin with a heuristic argument: It follows from 
the equipartition theorem that, in the limit fc^T hu 
where Q is the Debye frequency for the phonons of the 
vortex lattice, the lattice entropy is 



S = Nk B In 



k B T 



(5) 



where N is the number of phonon modes. For fixed to- 
tal number of 9 vortices, N v , piercing the system, the 
HQV lattice has twice as many phonon modes, N = 2N V , 
as the full quantum vortex lattice, since one full quan- 
tum vortex is topologically equivalent to a pair of HQVs. 
This means that the difference between the entropy of 
the HQV lattice and the full quantum vortex lattice is 
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As we are in the classical regime, we can expect the first 
term to dominate over the second term, so AS > 0. Con- 
cerning the energy, we can take an estimate from the 
'single particle' picture we had in the last paragraph: 



AE = N v k B T KT 



(7) 



where (nj p) c may be smaller or larger than 1 depending 
on the microscopic details (i.e. the ratio of the core ener- 
gies of full and half vortices) . From the thermodynamic 
expression T C AS = AE for a first order phase boundary 
and Eqs.([6| and ([7]), we find that the critical temperature 
for the transition from the full vortex lattice to the HQV 
lattice is 
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This equation indicates that for < k/ p — (k/ p) c <C 1, 
we have T c <c Tkt- 

In the next subsection, we will calculate the free energy 
of the phonon modes of the full and HQV lattices and 
justify our argument in this subsection. 



B. Phonon modes of an HQV lattice 

While it is possible to imagine circumstances in which 
complicated vortex lattice structures occur, we will here 
make the simplifying assumption that the HQV lattice 
consists of two interpenetrating triangular lattices, one 
of the (A6,Aa) = (n,ir) HQVs and the other of the 
(A8, Aa) = (ir, — ir) HQVs, translated by t with respect 
to each other. The full quantum vortex lattice in this 
picture is recovered in the limit r — > 0. In this paper, 
we consider the cases with highest symmetry, as shown in 
Fig. ([2]): a honeycomb structure analogous to graphene or 
an interlaced rhombus where a (A6, Aa) = (it, —tt) HQV 
lies at a midpoint between two (AO, Aa) = (tt, tt) HQVs. 
It is important to note here that, in an HQV lattice, 
there is a difference between 'intra-lattice' interaction U 
and the 'inter-lattice' interaction V: 



U = p 

2 M* H 

7T ft 2 

V = p 

2 M* H 



where i is the lattice point label for one of the two vortex 



lattices and r, 



t r 



Ti + r, and 



Ti'ji =v\— r'j = Yij. Note that whereas U is repulsive in 
both the 9 channel and the a channel, V is repulsive in 
the 9 channel but attractive in the a channel. 
The vortex equation of motion is 



■— = TThpZ X — Ti 
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(10) 



which is equivalent to the Hamilton equations of motion 
with Hamiltonian 



H v ({qi}) = U + V 



(11) 



where qf — \ptx~Kp and q\ = y/nhp ri tV are canonically 
conjugate variables. One important consequence of this 
is that the number of phonon modes for a vortex lattice 
is only half that of a usual crystalline system. 
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In order to obtain the phonon modes, we expand the 
equations of motion to first order in powers of the dis- 
placements, Sri, from the ground-state configuration. 
Because HQV lattice contains two vortices per unit cell, 
for each Bloch wave-vector in the first Brillouin zone, 
there are two phonon branches, an acoustic and an opti- 
cal. Near the Brillouin zone center, the acoustic modes 



involve the motion of the HQV pair center of mass 
fj = {t. l + r-)/2 while the optical modes are associated 
with the relative motion of HQV, Ari = — y[. (At 
large wave-vector, near the boundary of the Brillouin 
zone, these two motions are mixed, as shown in the Ap- 
pendix.) The long wave-length acoustic modes can be 
obtained from the equation of motion 
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and the corresponding optical modes from 
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This separation of relative and center of mass motion is 
justified in the limit where 
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Similar phonon dispersions were also f ound for the 
vortex-antivortex lattice in a 2D superflui d 33 * 34 *. 

Our result for dispersion of the acoustic branch justifies 
the arguments of Eq.(|6|. We show in Appendix A that 
for the acoustic branch, we obtain the same dispersion 
relation t hat Tk achenko (for neutral superfluid A — > oo)P^ 
and Fetter^SEZl obtained for the phonon modes of the full 



quantum vortex lattice: 



Tihn r 



Xdq 2 



2M* ^l + ^A) 2 ' 



(15) 



where n v = N v /Q (where fi is the total area of our vortex 
lattice) is the vortex density (where a HQV is counted as 
1/2 a vortex) and d is the distance between two closest 
HQVs of the same type. This derives from the fact that in 
the long- wavelength limit, the effective interaction for the 
acoustic branch is U+V. Because the acoustic dispersion 
is the same in both the full and HQV lattices, 
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where uj is the mean frequency of the optical branch. 
Note the similarity of this expression to the heuristic re- 
sult in Eq.([6]). 
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The optical branch dispersion we obtain from Eq.(13) 
justifies the modified version ofEq.(|| 



T r In 



kBT c 

hldr, 



T, 



KT 



(17) 



which we obtain from Eq.(f6|. This is because from the 
dispersion relation 
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we obtain \zbTkt ~ p/ n v 
implied condition for Eq.(17) - 



3> 1, which means the 
ksTc Hui - is satis- 
fied, unless k/p is fine-tuned to be very close to (/t/p) c : 
(k/p) — (k/p) c ~ n v /p. In order to avoid pathologies in 
the T — > limit, we treat the phonons as a collection 
of quantum harmonic oscillators, within the the Debye 
approximation for their spectrum, from which we obtain 
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which was used to plot Fig. ([lj; a here is a O(l) constant 
that can be obtained from the classical lattice energy 
calculation, which we will do in Section IV. 

However, Eq.(18) also implies that at T = the HQV 
lattice is dynamically unstable for any k > p. In the 
asymptotic limit that the vortex lattice is very dilute, 
where (k/p) c = 1, the T — transition between the two 
phases has the peculiar feature that it is first order in the 
sense that there is a discontinuity in the first-derivative of 
the ground-state energy at the critical point, but there 
is also a soft mode with a frequency that vanishes at 
criticality and no regime in which the HQV lattice is 
metastable for (k/p) > (k/ p) c . If corrections due to the 
core energies make (n/p) c < 1, the T = transition is 
more typical of first order transitions, in that both the 
full and HQV lattices have a regime of metastability be- 
yond the critical point. Below, we will show that ther- 
mal fluctuations can modify the phonon dispersion and 
increase the range of dynamical stability. Quantum fluc- 
tuations, which we have not considered seriously in the 
present paper, probably have a similar stabilizing effect, 
even at T = 0, and thus eliminating these pathologies. 



III. DYNAMIC STABILITY OF LATTICE 



A. "Self-consistent phonon" method for T > 



Hamiltonian which is quadratic in the deviations from 
the classical ground state configuration, 
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ij ab 



where Uo is the classical ground state energy, ij refer to 
lattice sites and ab are the coordinate components x and 
y. While at T — 0, K is equal to the appropriate matrix 
of second derivatives of H v , for T > we compute K 
according the the Feynman- Jensen variational principle 
so as to minimize the variational free energy, F': 



F'[{K?n] 
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where F is the exact free energy, (. . .} tr is the thermal 
average in the Gibbs ensemble defined by Ht r , and F tr 
is the free energy corresponding to H tr ^. With a little 
algebra, and with the use of Wick's theorem, it is easy 
to see that K is determined self-consistently from the 
relation 
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Such app roxim at ion has been used to calculate phonon 
dispersioi i 39 l 40 l. (We give a simple pedagogic example of 
this scheme applied to a single anharmonic oscillator in 
Appendix [B] ) 

To compute the first thermal corrections to the phonon 
dispersion, we expand U[{qf }] in powers of Sq^j, keeping 
terms up through fourth order. Then the self-consistency 
relation for the trial kernel can be expressed as: 
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where the derivatives of U are evaluated at <5<?£ = 0. 

In general, at low temperatures, since (5q^.5qf)t r °c T, 
the thermal corrections to the phonon frequencies are 
small at low T. The exception is as we approach a 
point of a zero temperature dynamical instability, where 
a mode goes soft, and the associated thermal fluctuations 
begin to diverge. The key point about the variational ap- 
proach (also illustrated by the simple pedagogic example 
in Appendix [B]) is that the fluctuations are computed 
self-consistently, so that the finite T regime of dynamical 
stability can even extend past the point of the T = 
instability. 



The dispersion relation obtained in the previous sub- 
section is modified at finite temperature by the anhar- 
monic part of the vortex interaction that we have so far 
ignored. In fact, we will show that when the anharmonic 
part is taken into consideration, the optical modes of the 
HQV lattice are stabilized even for k > p. 

To see this most simply, we treat the problem in a self- 
consistent phonon approximation. Here we define a trial 



B. Scheme for self-consistent dispersion formula 
for the HQV lattice 

The calculation of phonon modes proceeds much as 
before. Again, in the long wavelength limit, the acoustic 
mode is associated with the center-of-mass motion of the 
HQV pairs, while the optic mode is associated with the 
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relative motion. Following the same line of reasoning 
we followed in deriving Eq. 



12 and 13 the fluctuation 



operator can be expressed as a sum of contributions from 
the center-of-mass relative coordinate fluctuations in the 
long wavelength limit: 



(KM) 



(5Aq?8Aq%) (*A#*AgJ 



(SAqfSAql) , (5Aq?5Aq\ 



The effect of these fluctuations on the long wave-length 
acoustic modes is minor. However, on the optic modes, 
the fluctuation induced stiffening of the effective spring- 
constants makes the HQV lattice dynamically stable in 
the thermodynamic limit, at any finite temperature. In- 
cluding the corrections to the effective stiffness from 
Eq.(23l, the frequency of optical phonon at q = be- 
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where Arj is a particular combination of fluctuation 
terms computed self-consistently. It is straightforward 
but complicated to compute this quantity, as we show 
explicitly in Appendix [C] For instance, when the tem- 
perature is equal to the critical values at which the HQV 
vortex becomes thermodynamically stable, 
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Generally, the effect of the thermal fluctuations is largest 
anywhere that the T = phonons become soft, as 
here the thermal fluctuations are largest. Indeed, the 
impossibility of a dynamic instability at non-zero tem- 
perature can be proven by contradiction: were there a 
putative point at which the optic mode became gap- 
less, uj (q = 0) — > 0, the corresponding thermal fluc- 
tuations would diverge logarithmically with system size, 
(Ax? + Ay?) ~TlniV 

The same mechanism for attaining dynamic stability 
at non-zero temperature does apply in the neutral super- 
fluid (i.e. in the limit A — > oo). Possibly this signifies a 
dynamical instability to the formation of a HQV fluid, or 
it is possible that quantum effects or higher order non- 
linearities stabilize the HQV lattice. We leave this issue 
to be resolved in future work. In the next section, we will 
calculate the classical ground state energy for both full 
quantum vortex lattice and two HQV lattices of Fig. ([2]) . 
This will allow us to explicitly derive Eq.(7]) and also 
determine what structure the HQV lattice should have. 



IV. VORTEX LATTICE STRUCTURE 

We will first set up formulas for calculating the HQV 
lattice energy. We see from Eq.Q that the classical lat- 



tice energy per unit area is 
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for fixed magnetic field B = n„<I>o, where e' 1; e\ are the 
core energy of a single HQV and a single full quantum 
vortex respectively, 
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are dimensionless functions of the lattice structure, and 
(5e(r/£) vanishes for t/£ ^> 1, and accounts for the differ- 
ence between the core energies of the HQV and the full 
vortex when t/£ ~ 1. Note that Eo and E sp have im- 
plicit dependence on B since all distances between HQVs 
scale as roj, r^y cx y/~B. 

We find that the classical lattice Gibbs free energy, 
which determines the lattice structure, is proportional to 
the BKT temperature. This is because the quantity that 
can be tuned in a superconductor is the applied field 
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The thermodynamic relation between the interaction en- 
ergy and the Gibbs free energy shows that the latter 
should be proportional to the BKT temperature: 
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where B can be determined in terms of H from Eq. ( 28 ) 
and 



V A ) 



j2 ( r w_ 



\ A 



V A / 



K 



V A ) 

rof 
A 



(30) 



is another dimensionless function of the lattice structure. 

To determine the stability of different vortex lattice 
structures, we need to compute the Gibbs free energy as 
a function of the applied field H. For a given applied 
field, H , depending on the vortex lattice structure, there 
will be a different total magnetic field B. We can see 
this by inverting Eq.(28), assuming, as in Fig. [2j that 
the HQVs of the some type form triangular lattice. In 
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terms of H' cl = 8nei/$o, ignoring terms that are higher 
order in l/n„A 2 , 
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where, in terms of the exponential function (E±(x) = 
J x dte~ l /t), Jhqv is a r dependent function defined as 

Ihqv =E 1 (nnT 2 ) + Ei[™(r 0j - r 2 )] 
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7b (= 0.5772 • • • ) is Euler's constant, and 
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From Eqs.(29) and ( |3l[ ) (Si is evaluated in Appendix D), 
we obtain the Gibbs free energy 
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again ignoring terms that are higher order in l/7i„A 2 . 

We now can derive the vortex lattice phase diagram 
shown in Fig.Q by minimizing this expression for the 
Gibbs free energy. Eq. ( 34 ) shows us that n/p determines 
energy competition between the HQV lattices with differ- 
ent structures, as Ghqv is the only term in the equation 
that depends on the lattice structure. From Eq.(32), a 



straightforward numerical evaluation of the sums reveals 
that "fHQV smaller for the honeycomb structure than for 
the interlaced rhombus structure: jhqv ~ 0.4506 for the 
former and "/hqv ~ 0.5379 for the latter. Therefore, the 
honeycomb structure has lower energy for k < p and the 
interlaced rhombus structure has lower energy for n > p. 
Within our level of analysis, this tells us that if we ob- 
tain an HQV lattice at finite temperature through first 
order phase transition at k > p, the lattice will have the 
interlace rhombus structure. 

We can also determine the classical (k/p)c- For the 
triangular full vortex lattice, it is knowrF^ that 



Bfqv —H — H cl 
$ 



8ttA 2 



log 



4ir(H~H cl )\ 2 



IE - IFQV , 

(35) 



I 

and 
Gfqv/Q 



(H - H cl ) 2 
8tt 

32tt 2 A 2 

. 47T(g-g cl )A 2 

In v h 1-lE-lFQV 

(36) 

where H c \ = 47rei/<£>o- So the energy difference between 
the full vortex lattice and an HQV lattice is 



A(G/n) = {Ghqv ~ G FQV )/n 
t> (H-H cl 



In 



64tt 2 A 2 
4>„ 



K 

1 

P 



ir(H~H cl )e 
(H - H cl )AH cl 
An 



IE ~ IHQV + IFQV 

(37) 



where AH c i — H' cl — H c i . Note that since we are in the 
Hd < H < H c2 regime, $ » (H - H' cl )£ 2 . Thus, if 
we ignore the core energy, we have A(G/f2) < when 
p > k and A(G/f2) > when p < k. If we include the 
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core energy difference, the critical value will be 

1 - («/p)c 
f67rA 2 AiJ c i 



In 



n(H - H cl )£ 



2-lE-jHQV+lFQV 



(38) 



We see from this equation that (k/ p) c approaches f at 
lower vortex density, i.e. as H — > H c \. We also find, 



from Eqs(31|, (35), and (38), that there is a jump in B 



for any finite temperature transition to the HQV lattice: 
AB =B H Q V — B F qv 





K 




16ttA 2 


_P 


G)J 



In 



$ 



IE - 1HQV + 1FQV 



(39) 



Presumably, there are also quantum corrections to 
(k/p) c which we do not address here. 



V. DISCUSSION 

We can classify methods for detecting this entropy 
driven formation of HQV lattice into direct and indi- 
rect. The direct method would be vortex imaging, either 
through measuring the local magnetic field distribution 
using a scanning SQUID or Hall magnetometer, or with 
neutron scattering. One signature from vortex imaging 
that should be examined is the vortex lattice structure, 
since, as we see m Fig. p}, this is different for a HQV 
lattice compared to the full quantum vortex lattice. An 
indirect method would be to look for evidence of a first 
order vortex lattice to vortex lattice transition 
function of temperature. 

Throughout, we have assumed we have considered a 
system close to the zero temperature transition point, 
\n/p — (k/p) c \ <C 1, where the transitions take place at 
temperatures small compared to Tkt- At higher tem- 
peratures, other fluctuations will lead to a melting tran- 
sition^"^ f the vortex lattice itself. When the above 
inequality is well satisfied, these two sorts of transition 
can be treated separately, as they occur on very different 
temperature scales. Conversely, where it is not well sat- 
isfied, the physics of the full vortex lattice to HQV lattice 
transition may be altered significantly, and may even be 
preempted by melting of the full vortex lattice. 
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Appendix A: Phonon dispersion 



To calculate normal modes from Eqs.( 12 ) and ( 13 ), we 
set Ui = (C x ,C y ) exp[i(q-rj — iwi)] where u = f \ or £Arj. 
For both acoustic and optical modes, we can write down 
the normal mode equation of motion in the form 



a 



v 1 - 



»?(q) + C(q) a(q) 
"(q) »7(q) - C(q) 



which gives us the frequency 

^ 2 (q) = ^ 2 (q)-C 2 (q)-« 2 (q). 



c, 

(Al) 



(A2) 



As for these quantities rj, ( and a, we obtain for acoustic 
modes 

%(q) = YpliU) + (% r V)](l - e-^), 

Ca(q) = g^f Yy^o u ) cos + (%fV) cos 2<f> if ] 

x (l-e- iq - r *0, 
a (q) = — — y~][(fjjU) sui2<f>ij + (fiyV) sin 2^/] 



x (l-e- lq ' r ^), 
and for optical modes 



(A3) 



nhp 



x (l- e -' iq - r ^) 

+ -r- y^(Tij>V)cOS2(f> ir) 

nnp *■ — ' 

3 

a (q) = w^^2\{fijU)sin2<f) t] - (T ij ,V)sm2<j) ir } 



x (l-e^***) 



(A4) 
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where 



T 
t 



d 2 £P_ 
dx 2 dy 2 

d 2 



1 d 



dr \ dr 



1 d 



9 2 \ n ± n d 2 . n 1 d 
fTT^^T cos20 + 2— — -sin20 = r— . 
ox z oy z J oxoy dr \r dr 

(A5) 



and tan <f> = y/x. 

In the continuum limit, for the acoustic phonon pa- 



rameters we obtain 
hn v 



rdr / d(f)T 



and 

\( a (q) +ia a (q)\ = 



2M 

x (l_ e - i « rcos *) 

Trhn v 
= Af * A 2 . 

Trfet, (<?A) 2 
~ M* l + (gA) 2 



K, 



(D 



rdr[l - J (gr)]A 



2M * 7 
x (1 - e 
irhn v 
~ M*X 2 
nhn, 



rdr J d4>e m f 

iqr cos <fi \ 
/>oo 

\ rdrJ 2 (qr)K 2 



An 



M* 



(<A) 2 



1 + (<zA) 2 8 



2 i2 

—q a 



and for the optical phonon parameter 



(A6) 



(A7) 



hn v 
2M* 
tin, 



2M*X 2 



rdrd(j)(l-e- l ' lrcos4 ')T 
rdr I dcf>K ( 'J 



and 



|Co(g) + «a (g)| = 



hn v 
2M* 
hn v k 



-*iog : 



rdrdfa 1 ^^ - e~ iqrcos<t, )t 



log 



^((F 1 ) Tr/znt, k / 1 2 2 



, 1 - -q z d' 



2M* 



rdrd(j/T 



A 



2M* 



k l r 
-log 



A n 



log 



(A8) 



(A9) 



We obtain Eqs.( 15) and (181 by inserting Eqs.(A6), (A7|, 



and Eqs.(A8), (A9|, respectively, into Eq.(A2). To ob- 



tain finite temperature modification of dispersion, we 
therefore need to calculate the finite temperature cor- 
rection to these phonon parameters. 



Appendix B: Pedagogic example of the 
self-consistent phonon approximation 

We can show from an explicit example that at high 
enough temperature, the self-consistent harmonic ap- 
proximation around an unstable equilibrium may be more 
accurate than the approximation taken around a stable 
equilibrium. Our example is a particle in a two-well quar- 
tic potential, 

V, (Bl) 



with k, A > 0. At finite temperature, we can either take 
a self-consistent harmonic approximation around q = 



(which is the point of stable equilibrium for k < and of 
unstable equilibrium for k > 0) 

1 



H 



trl 



2m 



V 



1 Kq 2 , 
2 q ' 



(B2) 



where self-consistently K = -k/2 + ^/k 2 /A + 3\k B T/2, 
or (for k > 0) around a point of stable equilibrium 



Hi 



(B3) 



where self-consistently qo = ±yfc/A and K' — k + 
yjk 2 + 3AfcsT/2. The Feynman- Jensen free energy in 
these two cases are 

F[ = -k B T\nZ trl + (H- Utri)tri 



= — kf/T In 2itknT 



F' 



kgT In Z, 



tr2 



K 

(H — T-itr2)tr2 



= — kaT In 2Trkf>T 



knT- 



4A' 



(B4) 



respectively. 

Note that both solutions are at least metastable for 
all T > 0. For T > T', the harmonic approximation 
around q = has lower variational free energy than the 
one around q — qo, where T' is the solution of the implicit 
equation, 

k 2 . K' 2y/l + 3Xk B T'/2k 2 + 2 

= In — = In — . (Bo) 

2Xk B T' K ^/l + 6\k B T'/k 2 -I V ; 



Appendix C: Optical phonon stabilization at T > 

Here we examine the gap of the optical phonon (that 
is, the optical phonon frequency at q = 0) at finite tem- 
perature. We can expect the optical phonon stabilization 
to occur at the temperature where the optical phonon is 



gapless. Note that translational invariance ensures that 
acoustic phonon remains gapless at finite temperature. 
In order to calculate this, as indicated by Eq.(23), we 
apply the operator 



d 2 



1 3 2 



-i + i«<*> + <*">)(£ + £ 
+ !«*"> -<*■»(£-£ 



(SxSy) 



2 



dxdy 
(CI) 



on the vortex- vortex interactions U and V in Eq.(A4| 
We find 



2irhpr] (q = 0) = 



2irhp\{C,o + ia ){q = Q)\ 



2-K 2 h 2 pn v 
M* 



\Y^^ x l' + § yij') T if + i( 6x l' ~ s Vij>) cos 2^ + 2(Sx if Sy if ) sin 2^)77^7^ 



j 



2ir 2 h 2 pn v 
M* 



\{8Ax 2 + 5Ay 2 )Y,T 2 ,V + \ £>Aa^A % + SA yi SA yj )T 2 ,V 

3 3 

^ ^^((SAxiSAxj — SAyiSAyj) cos 2<j>iji + 2(5Axi5Ayj) sm2(/)ij')'flj>7ij>V 
j 

\Yl\^ x l + S V^3' + « 5 4 " s Vij) cos2^ tf / + 2<<y%/%jv) sin2^v)7;,v]7I J vV; 



K 2ir 2 h 2 pn v 
~p M* 



^1 E e2l0lJ ' [{5x 2 3 ,+5y 2 r )T ir + {(Sx 2 3 ,-Sy 2 3 ,) cos 20y, + 2{8x ij ,Sy ij .) ^2^}%']%^. 

3 

(C2) 



The thermal correction to |(£ + ia )(q = 0)| vanishes in 
the continuum limit because (Sxfj, + 8y 2 y) and (8x 2 j, — 
6y 2 3 -,} cos 2(j>iji + 2(5xiji5yiji) sva.2<t>iji do not depend on 
fay. Hence Eq.(24|. 

Among thermal corrections for 770(9 = 0), the biggest 
contribution comes from 



1 
4 

1 irh 2 pn v 
8 M*A 4 
1 irh 2 pn v 
8 M*A 4 



(^A^Ay 2 )^ 2 ^ 

(J Ax 2 + 5 Ay 2 ) / rrfr/v, 
(<5A* 2 + 5Ay 2 ) 



(C3) 



if the optical phonon is nearly gapless. Following 
Fetter, we find the thermal fluctuation due to the gap- 
less optical phonon with dispersion lo = rj d q/2 to be 



(SAx 2 + SAy 2 



2k B T r/ (q) 



4fc B T 



-Khpd1n v rj J q Trhpd 2 n v rj 



In AT,, 



(C4) 



As for other thermal fluctuation terms for 770(9 = 0), 
I Z) 3 -(<yAxi<yAxj + 5Ayi5Ayj)T?,V is a positive num- 
ber of comparable order whereas the rest of terms much 
smaller (~ lnn„A 2 ). Finally, we note that in calculating 
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the relative coordinate fluctuation, we have only com- 
puted optical phonon modes contribution, which is jus- 
tified as the In N v divergence originate from the infrared 
divergence of optical phonons. 



divergence and is convergent at short-wavelength. This 
second part is of order l/n v X 2 once we exclude the g = 
term. 

To calculate the first part, we first note 



Appendix D: Lattice energy summation 



Here, we evaluate Eqs. ( 27 ) and ( 30 ) . We will first com- 
pute ^ . K (r j'/X) through calculating^ 



(Dl) 



Using Fourier transform, we get 



= 2im v 

g 

= 2im v 



e i(g+k)-T 
(g + k) 2 + A-2 
o«(g+k)-x 



(g + k) 2 



2~Kn v X 2 



«(g+k)- 



(g + k)2[(g + k)2 + A-2]^ 



(D2) 



where g is the reciprocal lattice vector. Note that we 
have divided up the sum into the first part that is loga- 
rithmically divergent at both long- and short-wavelength 
and the second part that cancels out the long-wavelength 



I / y/4sirn^ poo 

(g + k) 2 =2 



ydye 



-y 2 (g+k) 



Jl/V47rn„ 
e -(g+k) 2 /47rn„ 



(g + k) 2 



2 I y dye- y (s+k) ". 

o 

(D3) 



Then we use Poisson summation 

47m„y 2 5>^ e ^ 2 (s+ k ) 2 = ^ e -*-(T- ri ) e -(T- ri )W 



3 



(D4) 

(this originates from J2 g e * s r = 2~2j <K r ' — r j)) t° derive 

(g+k)--r -y 2 (g+k) 2 



i 



g 

l ^ dy 



£ 



e ik.r 3e -(-r-r 3 ) 2 /V 



47m„ J y 



7TH .. ^ ^ 



m.. ^ — ' J 



(D5) 



Now we can evaluate "YljKoiroj' /X). The important 
step is to isolate out g = part: 



X)c ik - r ^o(A- 1 |r i -T|) 



2nn v > —r — 2im v X > 

^ (g + k) 2 ^ 



„i(g+k)-T 



(g + k)2[(g + k)2 + A-2] 



e j(g+k). 1 - e -(g+k) 2 /47 m „ I 

2nn v ^ (g + k)2 + 5 2^ e - t) ] - 2 ™- A £ 



i(g+k). 



g 

27m v e <k-T 
k 2 



(g + k)2[(g + k)2+A- 



-k /inn v _ 



A' 2 
k 2 + A" 2 



i(g+k)-r -(g+k) 2 /47rn„ i 

2™„ J] _ - 2 + - e^E^n^ - t)% (D6) 

g^O lg+ > ] 



where we dropped terms of order l/n v X 2 terms in the Combining this with Fetter's result^ 
last step. Taking k — ► 0, we obtain 



Y J K«{X-^-T\)=Y J K ( r ^- 



=27rn„A 



V A 

1 1 

2 + 2 



1HQV- (D7) 



= Trn^A 2 -^ ln47r?i u A 2 -i(l-7_E) + ^7 F Qy 

(D8) 
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gives us For the spin current part of the energy, we note 

V a =ATTn v \ 2 - l -\n^n v \ 2 - l -{2- lE ) e i(g+k)-x 

i 2 2 £ e .kr 3log( | r ._ T | ) = _ 2 ^^___ ) (DU) 

+ lflFQV + ^IHQV- (D9) J g 



We also find and obtain 

d 

^-U = 167rn„A" - 2, (D10) y 

2 """ t '^ 2 la ' 2'" rvr 2 

using dK (x)/dx = x[K (x) — K2{x)]/2. 



Z 1 =4n v \ 2 —, — r^S = 167rn„A 2 - 2, (D10) v 1 , C 2 1 , 1 1 / ni0 N 



1 M. M. Salomaa and G. E. Volovik, Phys. Rev. Lett. 55, 
1184 (1985). 

2 E. Babacv, Phys. Rev. Lett. 94, 137001 (2005). 

3 T. M. Rice and M. Sigrist, J. Phys.: Condens. Matter 7, 
L643 (1995). 

4 A. P. Mackenzie and Y. Maeno, Rev. Mod. Phys. 75, 657 
(2003). 

5 H.-Y. Kee, Y. B. Kim, and K. Maki, Phys. Rev. B 62, 
R9275 (2000). 

6 S. Das Sarma, C. Nayak, and S. Tewari, Phys. Rev. B 73, 
220502 (2006). 

7 S. B. Chung, H. Bluhm, and E.-A. Kim, Phys. Rev. Lett. 

99, 197002 (2007). 

8 J. Stenger, S. Inouye, D. M. Stamper-Kurn, H.-J. Miesner, 
A. P. Chikkatur, and W. Ketterle, Nature 396, 345 (1998). 

9 T.-L. Ho, Phys. Rev. Lett. 81, 742 (1998). 

10 T. Ohmi and K. Machida, J. Phys. Soc. Jpn. 67, 1822 
(1998). 

11 E. J. Mueller and T.-L. Ho, Phys. Rev. Lett. 88, 180403 
(2002). 

12 K. Kasamatsu, M. Tsubota, and M. Ueda, Phys. Rev. Lett. 
91, 150406 (2003). 

13 R. Barnett, S. Mukerjee, and J. E. Moore, Phys. Rev. Lett. 

100, 240405 (2008). 

14 P. Fulde and R. A. Ferrell, Phys. Rev. 135, A550 (1964). 

15 A. I. Larkin and Y. N. Ovchinnikov, Zh. Eksp. Teor. Fiz. 
47, 1136 (1964) [Sov. Phys. JETP 20, 762 (1965)]. 

16 E. Berg, E. Fradkin, E.-A. Kim, S. A. Kivelson, 
V. Oganesyan, J. M. Tranquada, and S.-C. Zhang, Phys. 
Rev. Lett. 99, 127003 (2007). 

17 E. Berg, E. Fradkin, and S. A. Kivelson, Phys. Rev. B 79, 
064515 (2009). 

18 E. Berg, E. Fradkin, and S. A. Kivelson, Nat. Phys. 5, 830 
(2009). 

19 Q. Li, M. Hiicker, G. D. Gu, A. M. Tsvelik, and J. M. Tran- 
quada, Phys. Rev. Lett. 99, 067001 (2007). 

20 J. F. Ding, X. Q. Xiang, Y. Q. Zhang, H. Liu, and X. G. Li, 
Phys. Rev. B 77, 214524 (2008). 

21 H. A. Radovan, N. A. Fortune, T. P. Murphy, S. T. Han- 
nahs, E. C. Palm, S. W. Tozer, and D. Hall, Nature 425, 
51 (2003). 



22 A. Bianchi, R. Movshovich, C. Capan, P. G. Pagliuso, and 
J. L. Sarrao, Phys. Rev. Lett. 91, 187004 (2003). 

23 D. F. Agterberg, Z. Zheng, and S. Mukherjee, Phys. Rev. 
Lett. 100, 017001 (2008). 

24 D. F. Agterberg and H. Tsunetsugu, Nat. Phys. 4, 639 
(2008). 

25 L. Radzihovsky and A. Vishwanath, Phys. Rev. Lett. 103, 
010404 (2009). 

26 N. Read and D. Green, Phys. Rev. B 61, 10267 (2000). 

27 D. A. Ivanov, Phys. Rev. Lett. 86, 268 (2001). 

28 A. Stern, F. von Oppen, and E. Mariani, Phys. Rev. B 70, 
205338 (2004). 

29 M. Stone and S. B. Chung, Phys. Rev. B 73, 014505 (2006). 

30 C. Wu, J. Hu and S.-C. Zhang, Int. J. Mod. Phys. B 24, 
311 (2009). 

31 J. Jang, D. G. Ferguson, V. Vakaryuk, R. Budakian, 
S. B. Chung, P. M. Goldbart, and Y. Maeno, Science 331, 
186 (2011). 

32 S. B. Chung, D. F. Agterberg, and E.-A. Kim, New J. 
Phys. 11, 085004 (2009). 

33 M. Gabay and A. Kapitulnik, Phys. Rev. Lett. 71, 2138 
(1993). 

34 S.-C. Zhang, Phys. Rev. Lett. 71, 2142 (1993). 

35 V. K. Tkachenko, Zh. Eksp. Teor. Fiz. 50, 1573 (1966) 
[JETP 23, 1049 (1966)]. 

36 A. L. Fetter, P. C. Hohenberg, and P. Pincus, Phys. Rev. 
147, 140 (1966). 

37 A. L. Fetter, Phys. Rev. 163, 390 (1967). 

38 R. P. Feynman, Phys. Rev. 97, 660 (1955). 

39 A. A. Kugler, Ann. Phys. (N. Y.) 53, 133 (1969). 

40 P. M. Platzman and H. Fukuyama, Phys. Rev. B 10, 3150 
(1974). 

41 A. L. Fetter, Phys. Rev. 147, 153 (1966). 

42 E. Frey, D. R. Nelson, and D. S. Fisher, Phys. Rev. B 49, 
9723 (1994). 

43 E. Babaev, A. Sudb0, and N. W. Ashcroft, Nature 431, 
666 (2004). 

44 E. Sm0grav, J. Smiseth, E. Babaev, and A. Sudb0, Phys. 
Rev. Lett. 94, 096401 (2005). 

45 A. L. Fetter, Phys. Rev. B 11, 2049 (1975). 



12 



